TELKOMNIKA Telecommunication, Computing, Electronics and Control 

Vol. 18, No. 5, October 2020, pp. 2708~2717 

ISSN: 1693-6930, accredited First Grade by Kemenristekdikti, Decree No: 21/E/KPT/2018 

DOI: 10.12928/TELKOMNIKA.v1815.5632 O 2708 


Accurate harmonic source identification using S-transform 


Mohd Hatta Jopri', Abdul Rahim Abdullah’, Rony Karim’, 
Srete Nikolovski‘, Tole Sutikno>, Mustafa Manap‘ 
“6Faculty of Electrical and Electronic Engineering Technology, Universiti Teknikal Malaysia Melaka, Malaysia 
Faculty of Electrical Engineering, Universiti Teknikal Malaysia Melaka, Malaysia 
Infineon Technologies AG, Germany 
‘Department of Power System, Faculty of Electrical Engineering, University of Osijek, Croatia 
`Deparment of Electrical Engineering, Universitas Ahmad Dahlan, Indonesia 


Article Info ABSTRACT 

Article history: This paper introduces the accurate identification of harmonic sources in 
the power distribution system using time-frequency distribution (TFD) 

Received Jan 15, 2020 analysis, which is S-transform. The S-transform is a very applicable method 

Revised Apr 16, 2020 to represent signals parameters in time-frequency representation (TFR) such 

Accepted May 1, 2020 as TFR impedance (Zrrr) and the main advantages of S-transform it 


can provide better frequency resolution for low frequency components and 
also offers better time resolution for high-frequency components. 
Keywords: The identification of multiple harmonic sources are based on the significant 
relationship of spectral impedances (Zs) that extracted from the Zzrr, consist 
of the fundamental impedance (Z7) and harmonic impedance (Zn). To verify 
the accuracy of the proposed method, MATLAB simulations carried out 


Harmonic source identification 
Power distribution system 


Spectral impedance several unique cases on IEEE 4-bus test feeder cases. It is proven that 
S-transform the proposed method is superior, with 100% correct identification of harmonic 
Time-frequency distribution source location. It is proven that the method is accurate, fast and cost-efficient 


to localize harmonic sources in the power distribution system. 


This is an open access article under the CC BY-SA license. 


Corresponding Author: 


Mohd Hatta Jopri, 

Faculty of Electrical and Electronic Engineering Technology, 
Universiti Teknikal Malaysia Melaka (UTeM), 

Hang Tuah Jaya St., 76100 Durian Tunggal, Melaka, Malaysia. 
Email: hatta@utem.edu.my 


1. INTRODUCTION 

As the power quality (PQ) issues caused mainly by harmonics directly affect the overall continuation 
of electric transmission and distribution network, the uninterrupted monitoring of electric power systems 
(EPSs) has become highly necessary for the utility industry [1, 2]. The harmonic distortion causes several 
adverse effects such as overheating of rotational machines, static machines, and current-carrying conductors. 
Furthermore, it causes premature breakdown or action of protecting appliances and harmonic resonance 
situations on the customers' EPS and metering imprecisions [2, 3]. 

In EPS, multiple harmonic sources can be found at the utility side, which known as upstream and 
downstream when harmonic sources located at the customer side with respect to a measurement point of point 
of common coupling (PCC) [4, 5]. Previously, numerous of harmonic source identification can be found in 
many kinds of literature. The earliest technique to localize the harmonic sources is utilizing power direction 
technique. However, this method used limited harmonic power-based indices, whereas mathematical analysis 
has shown the power direction method is not suitable for harmonic source detection [6, 7]. Furthermore, a state 
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estimation method is presented with the accomplishment of least square and Kirchhoff current law [8]. Yet, 
a study from [9], explained that the state estimation technique requires numerous measurement devices and 
unreasonable setup cost for the large power system. Thus, to overcome the previous constraint, the critical 
impedance method is introduced [10, 11]. The main disadvantage of this method is to know the value of source 
internal impedances of customer and utility, whereas it is difficult to obtain these parameters without switching 
tests [12, 13]. Consequently, a harmonic current vector method is introduced due to solving the previous 
method drawback. However, the measurement long-drawn-out and hard to analyze harmonic impedances at 
the customer side [14-16]. 

In the previous studies [17-20], the harmonic source estimation is discussed by applying the Bayesian 
estimation and harmonic state estimation (HSE). Unfortunately, this technique involves high computational 
complexity and the set-up cost of the distributed measurement system station is very expensive [5, 21, 22]. 
In real-time, accuracy, fast analysis and low cost are essential. Therefore, the limitation of previous analysis 
techniques can be solved by using signal processing techniques such as fast fourier transform (FFT), 
spectrogram and S-transform [2]. 

A harmonic source identification using FFT is introduced to provide fast analysis in determining the 
location of harmonic sources in the EPS. The significant relationship of spectral impedance, which are 
fundamental and harmonic impedance are used due to identify the harmonic location of the harmonic 
source [23, 24]. However, the drawback of the FFT technique, such as the signal information, for instance, 
amplitude frequencies and phases, cannot be acquired correctly. This is because of the leakage, aliasing effects, 
and picket fence formed by the FFT [25-27]. As a result of the FFT limitation, the spectrogram is presented to 
identify the harmonic source location [28]. Conversely, the spectrogram has two main disadvantages. First, its 
time-frequency resolution is limited by Heisenberg’s Uncertainty Principle. Secondly, both its energy location 
in the time-frequency representation and its resolution can vary with the window selected [29-32]. 

The S-transform is the advancement of wavelet transform (WT) and short-time fourier transform 
(STFT) by its variable window and phase correction, respectively. Similar to WT, it offers an improved time 
and frequency representation of a signal. Additionally, the scalable and movable Gaussian window properties 
of the ST can be utilized for superior recognition of any power quality events [2, 33-35]. Due to the excellent 
features of the S-transform in power quality analysis, it is a necessity to perform an analysis of identifying the 
location of the harmonic sources using this technique. In this paper, the S-transform is proposed to identify the 
location of the harmonic source in the distribution system with a single-point measurement approach at the point 
of common coupling (PCC). S-transform is a technique that uses time-frequency analysis to distinguish the signals 
in time and frequency representation (TFR) [36, 37]. The S-transform is proposed because it is able to provide 
provides better frequency resolution for low-frequency components and also offers better time resolution for high- 
frequency components. 


2. PROPOSED METHOD 

Harmonic current source type-load is selected as a frequent and widely use of harmonic producing 
loads in industrial [38, 39]. A time-frequency analysis which is used in this research is the S-transform. Three 
steps are proposed in this method. Firstly, measurement of voltage and current at the PCC. Secondly, these 
signals are analyzed using S-transform due to obtain the signal parameters in time-frequency representation 
(TFR) plane. Thirdly, the spectral impedance components are extracted from the TFR impedance and 
the significant relationship of these components are used in identifying the harmonic source location. 


2.1. S-transform 

S-transform combines the element of wavelet transform (WT) and the short-time fourier transform 
(STFT), which inherits the advantage of WT and STFT in signal processing. In particular, ST employs 
a moving and scalable localizing Gaussian window in the transformation process. The Gaussian window is 
as shown in (2). Mathematically ST can be defined as [29], 








ST (t, f) = E x(t) = oe 7 -ionft gt (1) 
g(t) = ere? (2) 
o(f) = + (3) 


If | 


where x(t) is the signal, t is the time, f is the frequency, g(t) is the scalable Gaussian window, and o({) is 
a parameter that controls the position of the Gaussian window. 
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2.2. Signal parameters 
The parameters of power quality signals are extracted from the TFR to provide the information of 
the signal in time. 


2.2.1. Instantaneous total non-harmonic distortion 
Total harmonic distortion, THD, is used to measure of how much harmonic content in 
a waveform [40]. The instantaneous total harmonic distortion of a waveform is defined mathematically as (4), 


ott Vh,rms (t)? 
THD(t) = L 


4 

Virms(t) ( ) 
where Vhr rms(t) is RMS harmonic voltage and Vr rms(t) is RMS fundamental voltage. The instantaneous total 
harmonic distortion of a waveform is defined mathematically as. 


2.2.2. Instantaneous total non-harmonic distortion 

Besides harmonic, a signal also contains interharmonic components that are not multiple integers of 
the power system frequency [27]. Referring to this as the instantaneous total nonharmonic distortion, 
(TnHD(t)), then for this parameter is written as (5). 


X | Vrms (t) =) Vh,rms(t)? 
TnHD (t) = = 


Virms(t) 


(5) 


2.3. Spectral impedance 
The spectral impedance (Zs) components are extracted from the TFR impedance (Zrrr). The Zrrr can be 
calculated using equation at each harmonic frequency and the values then plotted. The Zrrr is expressed as (6) [13]. 


_ Svt f) 
ZTFR = SEP) (6) 





where Sy(tf) is the TFR of voltage and S(t f) is the TFR of current. The spectral impedance consists of 
a fundamental impedance (Z1), whereas the harmonic impedance (Z+) is a harmonic impedance at the order 
number of harmonic. 


2.4. Implementation of proposed method 

In most recent studies, suggest the execution of the proposed technique can be realized as depicts in 
Figure 1 (a) and Figure 1 (b) using IEEE 4-bus test feeders, while the harmonic source is modeled as an 
inverter-based load, which is known as harmonic current source type-load [41-43], where N is a linear load, 
which is a resistor whereas H is a harmonic source. Besides, four specific cases are considered in this research, 
such as [15 , 44, 45]: 
- Case 1: No harmonic source in the power network system 
- Case 2: Harmonic source at downstream of the PCC 
- Case 3: Harmonic sources at upstream and downstream of the PCC 
- Case 4: Harmonic source at upstream of the PCC 






PCC PCC PCC 


l 
Downstream Case 1° Case 3 





Upstream 


(a) (b) 


Figure 1. (a) Upstream-downstream for case 1, (b) IEEE 4-bus test feeders for case 2, 3 and 4 


3. RESULTS AND ANALYSIS 
In this section, it is explained the results of the research and at the same time is given 
the comprehensive discussion. Results can be presented in figures and tables. 
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3.1. Case 1: No harmonic source in the power network system 

The linear load, which 1s resistive loads, are located upstream and downstream of the PCC. As can be 
seen in Figures 2 (a) and (b) show the measured voltage and current signals in the time domain. The peak 
voltage is 326.5 V. In contrast, the peak current is 66.5 A, respectively. Besides, Figure 2 (c) shows the THD(t) 
and TnHD(t) of voltage and both parameters are zero percent. Meanwhile, Figure 2 (d) presents the THD(t) 
and 7nHD(t) of current and both parameters are zero percent too. It shows that the voltage and current signal 
have neither harmonic nor interharmonic components. Figure 2 (e) shows the Zrrr of the PCC, while 
Figure 2 (f) presents the Zs. The most striking result from Zs, it is shown that the Z; at 50 Hz is of 4.8 ohm and 
no harmonic impedances exist. Thus, for case 1, the significant relationship between Z; and Z; at the condition 
of no harmonic source in the power system network, where for harmonic component, / 1s any positive integer 
whereas, for interharmonic, / is any positive non-integer. can be written as: 


Z, # Oohm (7) 


Zn # Oohm (8) 
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Figure 2. Case 1; (a) voltage signal in time domain, (b) current signal in time domain, 
(c) THD(t) and TnHD(t) of voltage, (d) THD(t) and TnHD(t) of current, 
(e) TFR of impedance (Ztrr), (f) spectral impedance component 


3.2. Case 2: Harmonic source at downstream of the PCC 

For case 2, the harmonic source is located downstream of the PCC. As can be seen in Figure 3 (a) and 
Figure 3 (b) show the measured voltage and current signals in the time domain. The peak voltage is 326.5 V. 
In contrast, the peak current is 66.5 A, respectively. In addition, Figure 3 (c) shows the THD(t) and TnHD(t) 
mean of voltage are 1.51% and 6.05%, respectively. In the meantime, Figure 3 (d) presents the THD(t) and 
TnHD(t) mean of current, which are 1.67% and 6.68%, respectively. The values of THD(t) and TnHD(t) show 
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that the voltage and current signals do have harmonic and interharmonic components. Figure 3 (e) shows 
the Zrrr of the PCC and the Zs is obtained by extracting the parameters from the Zrrr. Meanwhile, Figure 3 (£) 
presents the Zs that consists of the Z;. In contrast the harmonic impedance consists of Z275, Z375, Z600, Z700, Z825 
and Zooo, respectively. It is crucial to observe the relationship of Zs components due to identify the harmonic 
source location in the power system network. Interestingly, as shown in Table 1, the Z; value at 50 Hz, is 
always higher than any Z,. Thus, for case 2, the significant relationship between Z; and Z, at the condition of 
the harmonic source located at downstream of the PCC can be concluded and written as: 


Z, # Oohm (9) 
Zn < Zi (10) 


where for harmonic component, h is any positive integer whereas, for interharmonic, A is any positive 
non-integer. 
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Figure 3. Case 2; (a) voltage signal in time domain, (b) current signal in time domain, 
(c) THD(t) and TnHD(t) of voltage, (d) THD(t) and TnHD(t) of current, 
(e) TFR of impedance (Zrrr), (£) spectral impedance component 


3.3. Case 3: Harmonic sources at upstream and downstream of the PCC 
For case 3, the harmonic source is located upstream and downstream of the PCC. As can be seen in 
Figures 4 (a) and (b) show the measured voltage and current signals in the time domain. The peak voltage is 
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324.7 V, while the peak current is 44.6 A, respectively. Besides, Figure 4 (c) shows the THD/(t) and TnHD{(t) 
mean of voltage are 1.47% and 6.04%, respectively. In the meantime, Figure 4 (d) presents the THD(t) and 
TnHD(t) mean of current, which are 1.91% and 45.18%, respectively. The values of THD(t) and TnHD(t) show 
that the voltage and current signals do have harmonic and interharmonic components. Figure 4 (e) shows 
the Zrrr of the PCC, while Figure 4 (f) presents the Zs that been extracted from the Zrrr. Furthermore, 
Figure 4 (f) shows that the Z; is 5.36 ohm; in contrast the harmonic impedance consists of Z275, Z375, Z400, Z600, 
Z700 and Zooo, respectively. The characteristic of the spectral impedance, Zs are summarized in Table 2. What is 
surprising is that the value of Z; is always smaller than any harmonic impedances, Z,. It is apparent from 
the result that the Zs components can be characterized due to identify the location of the harmonic source. 
As a result, for case 3, the significant relationship between Z; and Z, at the condition of the harmonic source 
located downstream of the PCC can be written as: 


Z, + Oohm (11) 
Zy > Zi (12) 


where for harmonic component, / is any positive integer whereas, for interharmonic, A is any positive 
non-integer. 
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Figure 4 Case 3; (a) voltage signal in time domain, (b) current signal in time domain, 
(c) THD(t) and TnHD(t) of voltage, (d) THD(t) and TnHD(t) of current, 
(e) TFR of impedance (Zrrr), (£) spectral impedance component 


Accurate harmonic source identification using S-transform (Mohd Hatta Jopri) 


2714 O ISSN: 1693-6930 


Table 1. Magnitudes of the spectral impedance for Table 2. Magnitudes of the spectral impedance for 


case 2 using S-transform technique case 3 using S-transform technique 
Spectral Impedance Ohm Spectral Impedance Ohm 
ZI 5.36 ZI 6.5 
Z275 3.81 Z275 12.8 
Z375 3.68 Z375 11.8 
Z600 4.19 Z400 12.3 
Z700 4.34 Z600 12.2 
Z825 3.94 Z700 12.5 
Z900 3.68 Z900 12.9 


3.4. Case 4: Harmonic source at upstream of the PCC 

In case 4, the harmonic source is located upstream of the PCC. The voltage and current signals in 
the time domain are measured at the PCC are shown in Figures 5 (a) and (b). The peak voltage is 324.3 V, 
whereas the peak current is 30.65 A and the fundamental frequency is 50 Hz. These signals are used as 
the inputs for analyzing harmonic sources using the S-transform technique. It can be seen from the waveforms 
that the voltage and current waveforms are distorted because of the harmonic source presents in the network 
system. Furthermore, Figure 5 (c) shows the THD(t) and TnHD(t) mean of voltage, which are 1.47% 
and 6.04%, respectively. Moreover, Figure 5 (d) shows the THD(t) and TnHD(t) mean of current are 1.5% 
and 6.05%, respectively. The values of THD(t) and 7nHD(t) show that the voltage and current signals do have 
harmonic and interharmonic components. The Zs components are extracted from the Z7rr shown in 
Figure 5 (e). As shown in Figure 5 (f), the spectral impedance, Zs, consists of fundamental impedance, Z; at 50 
Hz, whereas the harmonic impedances such as Z275, Z300, Z375, Z400, Z600, Z700, Z900, respectively. 
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Figure 5. Case 4, (a) Voltage signal in time domain, (b) Current signal in time domain, 
(c) THD(t) and TnHD(t) of voltage, (d) THD(t) and TnHD(t) of current, 
(e) TFR of impedance (Zrrr), (£) Spectral impedance component 
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The Zs components relationship are the critical driving factor of harmonic source identification and 
the values of the components are summarized in Table 3. It can be seen that the Z; is always the same value as 
all harmonic impedances, Zn. Based on the results, it can be concluded that the Zs components can be used to 
identify the location of harmonic sources. Hence, for case 4, the significant relationship between Z; and Z, at 
the condition of the harmonic source located at downstream of the PCC can be written as: 


Z, + Oohm (13) 
as (14) 
where for harmonic component, h is any positive integer whereas, for interharmonic, h is any positive 


non-integer. 


Table 3. Magnitudes of the spectral impedance for case 4 using S-transform technique 


Spectral Impedance Ohm 
ZI 10.6 
Z275 10.6 
Z300 10.6 
Z375 10.6 
Z100 10.6 
Z600 10.6 
Z700 10.6 
Z900 10.6 


3.5. Performance of the proposed method 

This research provides an exciting opportunity to advance the understanding of harmonic source 
identification using the S-transform technique. The main contribution of this analysis is the significant 
relationship of the spectral impedance components, Zs such as fundamental impedance, Z; and harmonic 
impedance, Z; in identifying the location of the harmonic source in the network system. The spectral impedance 
characteristics are used to identify the location of harmonic sources and the performance of this proposed 
method are summarizes in Table 4. 


Table 4. Performance of harmonic source identification using spectral impedance components characteristics 
Case Upstream side Downstream side Spectral impedance, Zs characteristic % Correct identification 


1 N N Zı #0 ohm 100 
Z, = 0 ohm 

2 N H Z+ 0 ohm 100 
Zi > Zh 

3 H H Zı+# 0 ohm 100 
Zr > Z1 

4 N H Zı + 0 ohm 100 
Z= Zi 


4. CONCLUSION 

The main concern of this research is to identify the location of harmonic source that connected to 
the distribution system utilizing the S-transform technique. The main contribution of this research is 
the significant relationship of Zs components that obtained from the S-transform analysis in identifying 
the location of harmonic sources. Based on Table 4, the correctness of the proposed method is 100% for each 
case and the significant relationship of Zs for identification of harmonic source location are summarized 
as follow: 
- IfZ;40 and Z,= 0, there is no source in the network system. 
- IfZ,40 and Z; > Z;, there is a harmonic source located at downstream of the PCC. 
- IfZ;40 and Z> Z;, there are a harmonic sources located at upstream and downstream of the PCC. 
- IfZ;#0 and Z; -= Zn, there is a harmonic source located at upstream of the PCC. 
Consequently, considering the results, this method can be used accurately in distinguishing the location of 
the harmonic source. 
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